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Abstract: We study here the large N limit in the presence of magnetic monopoles in the 
Yang-Mills/Higgs model in three dimensions. The physics in the limit depends strongly 
on the distribution of eigenvalues of the Higgs field in the vacuum, and we propose a 
particular, nondegenerate configuration. It minimizes the free energy at the moment of 
symmetry breaking. Given this, the magnetic monopoles show a wide hierarchy of masses, 
and some are vanishing as 1/N. The dilute gas picture, then, provides an interesting 
framework for the large N analysis. 
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1 Introduction 

The subject of this paper is the analysis of some nonperturbative features of the Yang- 
Mills-Higgs (YMH) model in three dimensions for gauge group SU(N) and some impli- 
cations for the large N limit. 

Our YMH model is the direct generalization to SU(N) of the Georgi-Glashow model 
[|TJ . This was one of the first models providing a spontaneous breaking of a gauge sym- 
metry. 

The vacuum expectation value of the Higgs field provides a length scale and the size 
for magnetic monopoles, providing their stability. 

It was however argued in the known 't Hooft paper on oblique confinement pj that 
such configurations should exists also in pure Yang-Mills theory of strong interactions, 
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and play a major role in the mechanism of quark confinement with their condensation. 

Already the Georgi Glashow model in two dimensions shows analogous configurations, 
vortices ||, which account for the BCS theory of superconductivity 0. 

It was an idea from the early ages of QCD || |[| [Hj] that quark confinement could be 
explained with a picture analogous to that of type II superconductors but with the role 
of electric and magnetic charges interchanged (and for late issues like abelian dominance 



in confinement see HTl 



For nonabelian groups SU(N) the vortex solutions are not stable, mainly due to the 
triviality of the fundamental group, and this complicates the matter in that it requires 
the breaking to some subgroup with abelian factors. 

This is exactly realized in the YMH model, and the infrared behavior of gauge theories 



with compact gauge group [ 13| , was explicitly shown to hold in the context of the Coulomb 



gas approach to the ensemble of widely separated monopole solutions |L4[ 



This treatment is based on the fact that widely separated monopoles interact via a 
Coulomb interaction, despite the exact form of the energy is not known explicitly for 
generic parameters and also the space of moduli of the generic solutions has a nontrivial 
geometric structure |T5| . 

After these early years, lattice simulations have been largely employed to test the 
occurrence and condensation of monopoles in correspondence to the deconfinement tran- 



sition, in the pure Yang-Mills theory [16 



Because of the absence of the Higgs field, monopoles are of no fixed size, and thus 
some mechanism of generation of mass is necessary if monopoles are to be relevant for 
the confinement mechanism. This is indeed found on the lattice and there are arguments 
in the three dimensional continuum theory JT7|j . 

Let's recall also that the static four dimensional theory is equivalent to the YMH 
model with the limit of no potential. This is one of the main reasons to study the YMH 
model. 

Many interesting features are in the domain of the large N expansion, mainly: a) the 
classification of Feynman graphs according to the surfaces where they can be embedded 
gives also the order in 1/N, actually 1/N 2 , of the graph and this allows the contact 
with the string interpretation of gauge theories on one hand, and on the other with the 
random surface interpretation of zero dimensional matrix models; b) for suitable operators 
[ P^| ] the correlation functions "factorize", that is, are given by the disconnected part at 
leading order; c) finally N is not renormalized, in that it is a dimensionless external fixed 
parameter. 

Of course a special mention is due to the pure N = oo case, where one implicitly lets 
iV —>■ oo before removing any cutoff. This can bring the theory to have a different phase 
structure. There is a specific example in one dimension that is the Kazakov-Migdal phase 
transition. 

The second property above, factorization, shows that for the operators which follow 
it, the large N limit is a kind of semiclassical limit. Their fluctuations are suppressed, 
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and 1/N plays the role of h for them. 

One of the recent interesting phenomena, is the emergence of a new dynamic in the 
N = oo limit of some theories. These include affine Toda models, principal chiral fields 
in one and two dimensions, two dimensional QCD, not to mention the matrix models. 

In all these theories the effect of taking N large is to generate an infinity of states 
which coalesce to form collective excitations of some other, higher dimensional, theory. 

The emergence of a new dimension is the remnant of the matrix index of the diagonal 
fields. In these theories angular variables play an important role but are integrated, more 
or less explicitly, to leave an effective interaction for the eigenvalues. 



An example is the theory of the principal chiral field one the line|^0[, which is equiv- 
alent, for some boundary conditions, to two dimensional QCD with its interpretation in 
terms of two dimensional string. 

We come thus to our case of Coulomb gas of magnetic monopoles. It was constructed as 
a sum on the classical dilute configurations of monopoles, weighted with the determinant 
of gaussian fluctuations around them. 

With the Coulomb gas there is, via a Sine-Gordon transform, a dual representation 
and the possibility to achieve an estimate about the string tension for Wilson loops of 
large area. 

The string tension is related to the pseudo-mass of the monopoles by an exponential 
relation, which shows, like in the dual Landau- Ginzburg theory, that it is related to the 
density of monopoles. 

For SU(N) we have N(N — 1) species of minimal monopoles, with magnetic charges 
in different couples of U(l) sectors, and the coulomb gas can be generalized to this case 
[]21j] , as well as the Sine-Gordon transform. 

One has however to take into account the different species of monopoles that populate 
the model, and this was not done in [2lJ . Accordingly one has to consider the determinant 



of quantum fluctuations around the different monopole backgrounds. 
Up to this point the analysis is valid for any N. 

A possible new behavior comes from the large N limit, because there necessarily 
appears a distribution of masses which are present in the theory. 

These masses can be of order, a priori, in the interval from to N, and the physics 
at large N of course has to be very different from case to case. 

In this model the monopole pseudo-masses and the mass of the gauge bosons are 
governed by the differences of eigenvalues of the Higgs field at infinity, 0°°, so that all the 
model depends on its distribution of eigenvalues in the vacuum. 

In the standard picture of symmetry breaking 0°° can not be changed by any fluc- 
tuation, once the universe has formed. One simply fixes it. The modulus of 0°° gets 
renormalized and possibly shifted as with the arguments of the effective Higgs potential, 
but there is no indication on its direction in the Cartan space. 
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However in the course of our analysis we are lead to use the unitary gauge, because 
the physical degrees of freedom are explicit and the Higgs field is diagonal, and a curious 
phenomenon arises. 

The unitary gauge is somewhat singular, because its Faddeev Popov determinant is 
not defined in the continuum. This turns out to be the product at each point of the 
Vandermonde determinant constructed with the eigenvalues of 4>°°. 

This factor in the functional integral seems to provide measure zero for all configura- 
tions where some eigenvalues coincide, and to give a sort of repulsion of eigenvalues. 

This immediately faces with the problem that monopole configurations, which have 
necessarily points where the eigenvalues coincide, have all zero weight. 

Fortunately as the analysis is carried on, and still thinking that the theory is renor- 
malizable, we can show that once the quantum fluctuations of the massive gauge fields 
are taken into account, the Van-der- Monde ultralocal determinant is canceled almost 
completely. What remains is just the Van-der-Monde determinant of the eigenvalues of 

0°°! 

Still this term does not authorize us to think to a 'quantum lifting' of the degeneracy 
in Cartan directions, because again <f)°° is fixed at the "beginning of the universe" . 

However we think that at the epoch of its formation, the system is sensible to this 
term, and thus chooses the distribution of eigenvalues which maximizes it. 

We have found this distribution, which determines back the distribution of masses of 
gauge bosons and monopoles in the system. 

It shows a peculiar shape, and should bring peculiar consequences in the properties 
of the system. 

It is worth noting that the same distribution of eigenvalues of 'Higgs field' and gauge 
boson masses, is found |Z5 in the recent nonperturbative solution of SU(N) supersymmet- 



ric Yang-Mills theory in four dimensions (the M = 1 case), once the complete confinement 
is required. It is striking that the same distribution appears in the large N, and as its 
very origin is physically not known in that case. 

Coming back to our problem of the monopole gas, we have now a distribution to 
analyze the system and we could proceed. 

Unfortunately the very difficulty to the completion of the program is still there, be- 
cause the semiclassical sum needs some estimate with N of the determinant in the external 
field of a monopole. 



This problem is still unsolved despite many efforts |23| £4], £5|, |26[ and is nontrivial. An 
approach with Eguchi Kawai reduction which could provide at least the large N behavior 
is being studied. 

Let's summarize our path in this work. 

The first part of this work deals with some new large N ideas for the YMH model in 
three dimensions. Needless to say there is a large historical and scientific background and 
it is of course difficult to say something really new on these subjects. Nevertheless some 
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latest ideas on matrix models, For example, to our knoweledge, the study of Coulomb 
gases at large number of species is not investigated. 

After one finds a reliable method to estimate the functional determinant, it will be 
possible to draw definite conclusions on the large N monopole gas, which seems to promise 
interesting features, due to the competition of factors in the functional integral which takes 
place in the large N limit. 



2 Preliminaries 

We start by considering the Yang-Mills/HiggS3 model built on the gauge group SU(N), 
so it is a 3-dimensional euclidean theory of a gauge field = A a ^T a = A M • T and a 
matter field <p — <P a T a — cj> ■ T both living in the adjoint representation. 

Both are arranged in N x N matrix fields living in the algebra of SU(N). In all the 
following we will always denote with plain symbols such fields, like A^, <p, and with bold 
symbols the vector of their components along the algebra generators, like A M , <p above. 
Whenever the field will be diagonalized, its vector will have only the Cartan components. 

To fix the notation we take the normalization such that tx{T a T^) = \5 a b in the fun- 
damental representation. We will occasionally mention also the theory with T a in the 
adjoint representation, in next chapter, because they give rise to very different scenario 
of monopoles. 

The configuration space T is the space of functions from R 3 to the couple (A, <p) with 
finite action S. In this configuration space acts a continuous SU(N) gauge group: 

A^(x) -> oj~ 1 (x)d fl u(x) + goj~ 1 (x)A f ,(x)u(x) (2.1) 
(f)(x) -> w-\x)<i>(x)u){x) (2.2) 

which leaves the action invariant. 



Because we will focus on the large N expansion, it is necessary to adapt the parameters 
of the theory to this limit. 

It is the standard remark |27j that the large N limit is nontrivial only if the perturba- 
tive series remains finite, and this requires the coupling constants to be suppressed with 
powers of N. Rescaling the fields one can require all terms to be of the same order iV 2 . 

The action is then: 

r ( N N 1 

S = I I 2~2 tr ( G >' G W + yM^Z^) + V(4>) \ (2.3) 

with Vfj, = + g[Ap, •], = [2? M , V v \ = g(d^A u - dyA^) + g 2 [A^ A v \ and V(<f>) a scalar 
potential. 

This YMH theory in three dimensions can be seen as the static version of the four 
dimensional minkowsky YMH, or even of self-dual pure YM theory, at finite temperature, 
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where the A component assumes the part of the Higgs field (and V = 0). It has thus 
also direct phenomenological interest. 

So, to make contact with the four dimensional theories, one respects four dimensional 
renormalizability, and takes the potential to be a quartic polynomial in the traces* of the 
Higgs field: 

y(0) = A[tr(0 2 )- / x 2 ] 2 (2.4) 

From N power counting /i has to be of order N. 

Requiring the vacuum to have finite action (finite energy in four dimensions) V((f>) 
induces a spontaneous shift of the Higgs vacuum value from zero to some 0°°, and the 
gauge group G breaks down spontaneously to a subgroup H which leaves 0°° invariant. 
This is called the little group* . 

T is divided in disjoint sets classified by the winding of the two-sphere at infinity into 
the coset G/H: elements of ^(G/H) which is isomorphic to tt\(H) because ivi(G) is 
trivial. 

In case the vacuum Higgs field has all different eigenvalues the little group H is the 
maximal abelian subgroup [/(l)^" 1 , and the classification has N — 1 topological quantum 
numbers: Z^" 1 . Moreover thanks to the vanishing of tt 2 (G), this classification is gauge 
invariant, because gauge transformations of (p^ are homotopic to the identity. 

Then, in the spirit of the semiclassical approach, one considers the minimum of the 
action in each set (j)(x), and expands <\> around it. One has a pointwise breaking of the 
local gauge group G down to the little group which leaves <j)(x) invariant. All gauge fields 
not belonging to H acquire a mass, while those in the little group remain massless. At 
points where the Higgs field has all different eigenvalues, the breaking is maximal and the 
little group is U(1) N ~ 1 . 

As 't Hooft shows ||, the manifold of points where two eigenvalues of 4> coincide has 
dimension d — 3, that is, in three dimensions, consists of isolated points. There the field 
configurations have the properties of magnetic monopoles, and in the next section we will 
review also, in the case of SU(N), how the topological number which classifies the Higgs 
field represents the magnetic charge of the configuration under the broken symmetry 
group H. (modulo equivalence under the Weyl discrete symmetry). 

The other terms in the action, namely the pure gauge and the Higgs kinetic terms, 
pose no obstruction*, and configurations of arbitrary winding can be explicitly constructed 

MM- 

*The choice of a function which is symmetric in the different Cartan directions is not the only one 
f2Sj but is the more natural. 

'''The little group H is always of the form T" x G', where T" is some abelian group and G is a simple 
subgroup. In any case it always includes a U(l) subgroup, called the electromagnetic group and generated 
by the Higgs field itself. 

*This happens in the gaugeless limit, where, turning off the coupling to the gauge fields, the Higgs 
kinetic term forces <f> to have the same direction at infinity, thus giving winding zero. 
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Before reviewing in the next section the classification of classical solutions for SU(N), 
we make a remark that we will need later valid for all the configurations of nontrivial 
winding. 

Following <p smoothly in all the space, one necessarily meets points where it has at 
least two coinciding eigenvalues, because otherwise the winding would have disappeared. 

This statement is obviously gauge invariant so it is true even in non regular gauges 
like the unitary gauge. 



3 Classical monopoles 



In this section we review the classification of monopoles [|, [31], [32], |33|, [34], |35], |3G|, [37|] that 



arise in 3 dimensions for the YMH model with group SU(N), that we will need later. 
The action is: 

1 ._ „ A 



S = N / d 6 x 



^tr(G^G^) + |tr(P M 02V>) + ^(tr(0 2 ) - /i 2 ) 2 



(3.1) 



First of all one introduces what are called point monopoles, singular configurations 
satisfying the equations of YMH plus the requirement of minimum action: 

V li G (tt , = ^) = V^ a = (3.2) 

The last two equations impose that one gets everywhere zero contribution to the action 
from the Higgs sector, so that we are in what is called Higgs vacuum. Some gauge fields 
have to be zero and the others allow for nontrivial solutions localized near isolated points. 

To see this it is useful to choose the abelian "unitary" gauge (ft G TC (with TC the 
Cartan subalgebra). We will write = <fi- T, where T is the vector of (iV — 1) commuting 
generators of TC in some faithful representation and is a vector in R^ -1 . 

In this gauge the above equations imply that the Higgs field is constant in all space 
4>(x) = (p^ and that has nonzero components only in the algebra of the little group of 
0oo, i.e. only if [A, </>«,] = 0. 

For 0oo with generic eigenvalues (all different) the little group is simply [/(l)^ -1 , so 
that only the abelian, Cartan, gauge fields survive in the vacuum. 

Here 0oo induces perturbatively a mass term for each gauge field that is charged 
with respect to it. For each root ct of SU(N) the mass of the charged A^ fields is 
m w a = 914*00 ' a \- At the same time the Cartan gauge fields decouple completely from 
the Higgs fields. 

As V{4>) is flat in all gauge directions, the V'(<p) = constraint only fixes the modulus 
of the vacuum Higgs field: tr^ = /i 2 . 

We will assume at this point a vacuum 0^ with all different eigenvalues. This is 
a gauge invariant statement and we will see that this vacuum configurations should be 
preferred by the system itself, in the large N limit. It could also be imposed by some 
gauge invariant external source. 

(checked to here) 
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3.1 Point monopoles 



The G^y field equation now allows for abelian U(1) N 1 solutions with a singular Dirac 
string [Q: 

= T • q/^ = (1 - ^/r^tan^/x) (3.3) 

= <7 1 ■ q | | 3 q e k 

q is the nonabelian charge of the monopole, but as we are in the unitary gauge, it belongs 
to XJ{X) N ~ X : it is for now an arbitrary vector in R 7 ^ 1 . 
As z — > — oo we have the asymptotic form: 

— 2T • qd M $ (0 < $ = tan _1 (y/x) < 2tt). (3.4) 

Observing the phase of a loop around the string we get a realization of 7Ti(f/(l) iV_1 ) 
and obtain the admissible monopoles: the generalized Dirac condition 

e 4*iq.T = ^ ( 35 ) 

This condition restricts the possible charges q to belong to a lattice in R^ -1 , in fact 
q has to be reciprocal to each weight of the representation chosen for the T: for every 
weight nij, 

q ■ mi = y n G Z. (3.6) 

The lattice of charges depends thus on the representation chosen for the T's, calling 
in the game also the global properties of the representation of the gauge group. 

It can be more or less dense depending on the modulus of the highest weight of the 
representation. 

One now introduces the (dual) co-roots, at* = oc-jcxi ■ cti, (where cti are the simple 
roots). For each weight they satisfy the relation m- a* = n/2, so that they are reciprocal 



to the weight lattice. The coroot system defines what f3T|] have called the dual group. 
For SU(N) the dual group is isomorphic to it, denoted SU*(N). Moreover for roots 
normalized to unity the coroot lattice coincides with the root lattice. 

An immediate consequence of the relation with the weight lattice is that the coroots 
(and also their multiples) are always between the possible magnetic charges q. Monopoles 
in the adjoint representation of the dual group are thus always present. 

Usually one studies the simplest cases of the fundamental and adjoint representations, 
for the gauge field. 

For T in the fundamental representation (T = |A) the weights are the fundamental 
ones and the magnetic lattice is just the coroot lattice 

q = J2 n i a i n { eZ (3.7) 

i=l...jV-l 



8 



The monopoles of minimum charge transform in this adjoint representation. 

The picture represent the coroot lattice for SU*(3), its generators as black circles. 
The small triangles represent the fundamental monopoles that arise for adjoint gauge 
generators^. There is thus a nice duality: for gauge variables with fundamental (adjoint) 
generators, the minimum monopoles transform in the adjoint (fundamental) representa- 
tion of the dual group. 




• Minimum adjoint charges 
<> Minimum fundamental charges 



A remark is due for the Weyl group which acts in the weight lattice and thus also 
on the magnetic charges. It is generated by the reflection with respect to the planes 
orthogonal to the roots, and sends every lattice that we have considered into itself. The 
action can be seen as a reflection also on the magnetic charges. 

Seen on the Cartan generators, it simply exchanges the diagonal entries (as can be 

/ \ 

easily seen in the (overcomplete) basis 2(Tj,)^ = (Sik$ji — SuSjk) — )• 

I 'J 

On the fields, any Weyl action is equivalent to a global gauge rotation with respect 
to the generator E a + E- a , where a is the relative root, and this shows that monopole 
configurations related by Weyl symmetry are gauge equivalent. 

This has implications on the type of monopoles: for the usual case of gauge variables 
in the fundamental representation the minimum charge monopoles are all related by Weyl 
t r ansf ormat ions . 

In the adjoint representation instead monopoles are classified by the dual fundamental 
weights which are not all Weyl-equivalent: they are divided in classes according to the 

§In fact an other case mentioned in the literature j2l[] is that of gauge variables with generators in 
the adjoint representation. The nonzero weights are in this case the simple roots, so that the reciprocal 
lattice coincides with the weight lattice of the dual group (m* • a = n/2): 

q = ^ n»m* n, £ Z. (3.8) 

i=l...JV— 1 

One can say that the minimum charge monopoles transform now in the fundamental representation of 
the dual group. They are shown in the picture as small triangles. The monopoles that arise in this case 
include also the previous adjoint charges as combinations of minimal monopoles, although the lattice 
generated by them is not shown. 
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N — 1 (nontrivial) elements of 7Ti(SU(N)/Zn). Weyl reflections only act within these 
classes. 

An example for the dual SU(3) is shown in the picture, the fundamental co-weights are 
represented by the the small triangles, the right pointing giving the [3] representation, the 
left the [3*] . The Weyl transformations are the reflections with respect to the long-dashed 
lines. 

One sees that the fundamental monopoles come in two triplets invariant under Weyl, 
while the usual adjoint monopoles come in a whole sextet. 



3.2 Regular monopoles 

One notices first from (|3T5[) that if the C/ JV_1 (1) path e 2q ' T * (0 < $ < 2ir) can be con- 
tinuously gauged away in SU(N) to the identity, then the Dirac string will decrease of 
intensity and disappear. The gauge transformation needed to do that is necessarily non- 
costant, so that one will end up with a noncostant Higgs field. If this can be done, the 
point monopole is the basis for a regular one with finite energy. It is clear that there is a 
great freedom to construct these regular solutions of the equations of motion. 



The spherically symmetric solutions have been studied in detail in |32[ with a classifi- 
cation of all the possible magnetic charges from which one can determine a finite action 
configuration. 

The charges which admit spherical solutions are given by q = q' — q" where q' and 
q" are the roots of two embeddings of SU(2) in SU(N) and q" must also be in the little 
group. One SU(2) q' is needed to rotate the Higgs field to a radial gauge, and the other 
is a remaining freedom to define the spherical gauge configuration. 



Since the factor in (|3.5| ) is a loop in the chosen representation of SU(N), the process is 



clearly possible for any charge if and only if the generators are in a faithful representation 
like the fundamental one^. 



The Weyl group relates monopoles of different charges by a global gauge transfor- 
mation. This is a difference with the abelian case where different charges classify gauge 
inequivalent monopoles. This peculiarity of non-abelian theories follows mainly from the 
simple-connectedness of SU(N), but also from the fact that colored flux lines are not 
gauge invariant (but covariant) and can thus be deformed or changed of color by gauge 
transformation. 

^In the case for example of the adjoint representation, there are ./V inequivalent loops which join the 
identity to the elements of the center of SU(N). So, for N — 1 kind of point monopoles, the string is 
impossible to remove and they are genuine Dirac monopoles. An other accident of this case is that the 
minimum charge monopoles are in this case N(N-l) (the weights of all fundamental representations), 
while the nontrivial elements of the center Zjv are N-l. The minimum charges can thus be divided in 
(N-l) sets of N elements, and the Weyl group acts only within each set. 

All this fundamental monopoles, associated to the nontrivial paths in SU(N) around the Dirac string, 
cannot be made regular. 



10 



Another important notice is that the charge quantization condition (|3.5| ) does not 
depend on the vacuum Higgs field 0^, and this is an advantage of the unitary gauge. 

Regular monopoles instead do depend on the 0oo boundary condition, because they 
are constructed transforming the Dirac string into a varying Higgs field. 

As we are going to assume 0oo to have all different eigenvalues, the charge q", belonging 
to the little group, can only be zero, so that the magnetic charges of spherically symmetric 
monopoles will coincide with q roots of SU(2) embeddings" . The figure shows them for 
SU(3), for iV > 3 the pattern is much more complicated. 

Possible spherical monopoles for SU(3) 

• 3 -> 2+1 embedding Q=(l/2,-l/2, 0) 

O 3 -> 3 embedding Q=(l, 0, -1) 

□ 3 -> 3 only with degenerate <|> Q=( 1 ,-1/2,- 1/2) 




Only the first two cases are possible for nondegenerate 0oo. For this cases the solution 
are well-known 't Hooft-Polyakov monopoles: 



A*(r,q) 
0(r,q) 



TqD^ + K{m w r) [E^iie + + E^{-i§ + 
^ • (T - qT q ) + if(m H r)(0 oo ■ q)T q 



(3.9) 



T q = q • T with some £" q and £"_q, give the SU(2) subalgebra. 

Changing to a regular gauge one can "smear out" the string in and leave an isolated 
singularity at the origin. 



3.2.1 Masses 

In the last expressions, mw is the mass of the two charged gauge bosons in the chosen 
sector of the monopole, function of the Higgs vacuum: my/ = g\<p°° ■ q| = — (pfl; 

m H is instead the mass of the Higgs field "in the sector of T g ": m\ = 2-^(0°° — 0°°) 2 . 

They regulate the exponential decay of the massive field components of the monopole 
as well as its total classical action. We have in fact 

Sd = m ^ Ci ^ } . ( 3. 10 ) 



"Monopoles of higher charge, empty sites in the picture, are not realizable as single three dimensional 
spherical configurations, although the topological argument indicates the existence of some extrema of the 
action. Some can be constructed as multimonopole-like configurations which possess discrete symmetry 
groups, with stability given by the Higgs attraction. Such configurations have been found to exist with 
tetrahedral, octahedral (but not icosahedral symmetry, for a late reference, see ]39| and therein). 
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The function C is found [EE] to approach the value Air when the argument goes to 
zero, which for 577(2) is called BPS limit. 

Considering the N(N — 1) unit charge monopoles, we see that, fixed </>°°, their pseudo- 
mass ranges in the interval < S c i < AnNfi/g. 

We see that most of the properties of the various objects are ruled by the m\y of the 
gauge bosons in the relative SU (2) sector, so that all depends on the set of Higgs vacuum 
eigenvalues (pf. 

The constraint \<j>°°\ — [i — 0(N 1 ^ 2 ) has important consequences on the set of masses 
that are present in the model in the large N limit, because, as we will see, necessarily 
there are masses that become small at least as 0(1/ N). We will also find masses of order 
0(1/N 2 ). 

All the monopoles of charge greater than 1, have higher mass proportional to their 
charge and they are expected to dissociate into smaller constituents. Hence their contri- 
bution to the infrared region is negligible and one can safely discard them. 

We will instead use approximate solutions built by superimposing an arbitrary number 
of minimal monopoles at large distances. 

They are constructed easily in the unitary gauge and are then regularized by means 
of a procedure similar to the one for the single monopole. We just need the proof of 
existence of the gauge transformation needed to do that because we will work in the 
unitary gauge. 



Ap(x) = T l n a =iA tM (x-x a ,q a ) 

{n) (^)=Ea=l(0(^-^,qa)-0 OO )+0° O { ' ) 

They depend only on the parameters of the n single monopoles {x a , q a } and are good 
solution to the equations of motion for distances much bigger than the monopole sizes. 
This approximation is called dilute gas. 

The interaction of monopoles simplifies drastically for large distances compared to the 
monopole size and remains function of the relative distances only. 

In fact the action for such dilute configurations is found to be approximated by the 
self-action of each monopole plus a Coulomb-like monopole-monopole interaction term 

ph. 



4 Quantum fluctuations 

Coming to the much more ambitious program of quantizing the theory, the only feasi- 
ble way is for now the semiclassical expansion, leaving apart the S-duality which solves 
the supersymmetric theories, or the first order (BF) formalism, which also falls back to 
semiclassical quantization. 
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If semiclassical it is, it should show the correct result. 

However even this treatment at one loop order is a nontrivial task, because involves 
the calculation of functional determinants in classical backgrounds. 



There are some simplifications in the BPS limit |38|, [JOf of supersymmetric theories, 
because the three-dimensional configurations represent four dimensional selfdual back- 
grounds and there are useful relations between fermion and boson determinants in this 
case, but we remind that the physics in this limit is drastically different. The main reason 
for this is that the Higgs field is massless and thus gives a further long range interaction. 
We thus want to consider A > 0. 

A major progress was made by Polyakov. In analogy with other simpler models he 
applied the semiclassical quantization to the 577(2) monopole solutions and carried 



the program to evaluate the Wilson loop resumming the semiclassical expansion. 

His treatment, for the compact QED, shows that the area law for the Wilson loop 
emerges from this semiclassical treatment precisely because of the condensation of mono- 
poles. 

Later in |2TJ Das and Wadia have reached the same conclusion for the problem with 
577(3) gauge group. Recently also for the case of pure YM 3 they have argued that 
confinement arises from the gas of monopoles, using nonperturbative results from the 
theory of the three dimensional Coulomb gas |17 |. 



So we have, at our disposal, the minima of the classical action that are strongly 
supposed to give the main contribution to the functional integral, and taking into account 
the gaussian fluctuations around them, one can introduce a semiclassical sum. 

We will see that important nonperturbative features of the model are reproduced by 
this approach. 

4.1 Semiclassical program 



In [14], |2TJ] the semiclassical quantization of a system with monopoles is approached 
through a grand canonical ensemble of magnetic particles. The sum on all configurations 
gives the partition function in a nonperturbative way, and after a generalized Poisson 
transform, a saddle point technique can be applied. 

The gauge fields should be integrated perturbatively at one loop, taking into account 
the regular monopole backgrounds as nontrivial minima of the action. 

In this approach many approximations are to be taken carefullly: 

First, the gas of monopoles is assumed dilute, thus considering just the coulomb part 

of the monopole-monopole interaction and simplifying the classical and one-loop effective 

action. 

After the generalized Poisson& Sine-Gordon transform has been used, usually one is 
limited to the minimum charge monopoles, assuming the higher ones dissociate rapidly**. 



**One could try to consider also the higher charges. From the Wilkinson Goldhaber analysis, that 
spherical monopoles have limited charge (for limited N, |q| < N — 1 ). Moreover their mass grows 
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4.2 Grand canonical ensemble 



The partition function of pure YM in 2+1 dimensions is transformed into the sum on all 
configurations made of any number of monopoles of charges {q a } and locations {x a }. For 
SU(2) it is: 

oo tn 



/ -27r Qalb 
IJ d 3 x a e^ ^ {q a = ±1, ±2, . . .) 



{qi-q™} 



O 



where £ is the classical and one loop contribution to the action of the one monopole 
configuration [fL4jl . 

For SU(N) this expression is no longer valid because the weight of each monopole 
depends on the charge, £ = £(q). Hence one has: 

OO 1 

3 = E zyGn (4-1) 

n=0 /t- 



d z„ C(q a ) • e~^~ i^-^i (q a in £/te root lattice) (4.2) 

iqi-qn) a=1 



The partition function Z can be nicely reexpressed as a functional integral using the 



Sine-Gordon transform 41 



Z = J V X e-^I d3x [^ 2 -^i M2 ^ eicii - X }. (4.3) 

The mass M 2 comes from the weight £: M 2 = 32ir 2 C,/g 2 . x is a iV — 1 components scalar 
field whose propagator is just the coulomb potential, and the sum on i is the sum on 
all the possible magnetic charges of one monopole (the magnetic lattice). The last term 
in the integral is the functional generator of the multi-charge configurations. For the 
symmetry of the minimum charges, the roots q^ = — q_j, the last term becomes a cosine. 

This representation of the coulomb gas has to be understood in a perturbative sense, 
because it is just the perturbative expansion which reproduces, diagram by diagram, the 
dilute gas. In this spirit the x field configurations have to vanish at infinity, although 
there is formally an infinite 'zero-point' energy of the vacuum x = 0, which disappears 
with the normalization. 



linearly with the charge while their size gets linearly shrinked. 
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4.3 Wilson loop 

One can then succeed to evaluate gauge invariant operators like the Wilson loop Wc- 



W, 



c 



N 



(4.4) 



for any contour C in three dimensional space-time. 

If we take into account the form of in the unitary gauge ( |3.cS| ) for each monopole 
of charge q a and location x a , (and after using the Stokes theorem) we can rewrite the 
Wilson loop as an external source for the configuration of charges {q a }: 



W, 



c 



N\ 

— /Tre^'Sa^^) 
N \ 



T)(x) 



d 2 a,(y) 



(g - y) f 

\x — y\ 3 



(4.5) 



Instead of the flux through the loop of 
the magnetic field produced by the q a charges, 
one thinks to a potential (tj(x), in the pic- 
ture) which acts on the charges q a at points 
x a , produced by a dipole layer on the sur- 
face S spanning the loop^. The problem 
of evaluating the Wilson loop in the func- 
tional integral is reduced to the average of 
the canonical ensemble under the action of this external potential. 

In the limit of very large loop this potential is constant on the two sides leaving just 
the discontinuity, which is of crucial importance. 

From the color trace the potential takes each direction along the fundamental weights 



T1(X) 






2n 










X 




-27T 





m. 



W, 



c 



(4.6) 



It is straightforward to evaluate this operator in the coulomb gas ensemble ( |4.1| ), and 
the result is, after a shift in the x field: 



W, 



c 



(4.7) 



4.4 Saddle point solution 

The nongaussian integral above is solvable by saddle point, remembering that M 2 is 
asymptotically small. 

The saddle point is given by the equations, one for each m a : 

-d 2 Xa + m a d 2 V = - £ M 2 (q,) qi sin(q,, ■ Xa ) (4-8) 

half i 

TTThe dipole density is unitary and in direction of d 2 ^ orthogonal to the surface. The potential is in 
practice the solid angle of the loop seen by the charge at x a . 
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(The sum on % runs now only on one half of the magnetic lattice because the negative 
symmetric part is included in the sine). 

The term m a d 2 r] imposes a discontinuity of the solution in 'internal' direction m a 
and together with the constraint of perturbativeness x(±oo) = 0, this leaves only the 
solutions of constant direction x( x ) — m aXa( x ) (no sum). 

Taking the product with m a , one gets the scalar equations 

(d 2 x a - d 2 7]) -— 1 = M 2 (qi)m a ■ q { sin(q, • m a x a ) (4.9) 

half i 

(where we used the fact that for the weights of the fundamental repr. m 2 = 4^-). 

At this point one has to specify which magnetic charges q» to use, and we see from 
the Sine-Gordon representation that higher charge monopoles give rise to a generalized 
Sine-Gordon potential which is of lower magnitude, even if it has shorter periods. This 
means that higher charges only perturbate the potential given by the minimal charges. 

We limit then the sum over i to the minimum magnetic charges, i.e. the co-roots of 
the first picture. 

We have to evaluate the scalar product m a • q$ of a fundamental weight with the 
adjoint weights. 

To this aim we remember that the roots qj are N(N — 1), and that given a fundamental 
weight m a , (N — 1)(N — 2) of them are orthogonal to it, while N — 1 have scalar product 
1/2 and the others (N — 1) (negative symmetric) have scalar product —1/2. 

It is then sufficient to limit the sum to the (N — 1) cases all giving result 1/2 for the 
scalar product: 

d\ a = d 2 rj + sin( Xa /2)(iV - l)Ml (4.10) 

where we have introduced the averaged M 2 given by 

«=(^T) MM.)- (4.11) 

It is this quantity which carries information on the physics that we obtain in the large 
N limit. Each M 2 (q) is, in Coulomb gas language, the fugacity of the monopole species 
q. M 2 is then directly related to the average density of monopoles, which is strongly 
believed to be the order parameter for confinement. 

We note also that many monopoles are mutually neutral and that M 2 is the average 
in the (N — 1) sectors that are not orthogonal to m a . To evaluate the Wilson loop we 
should also average on m a , finally. 

The solution of ( |4.10| ) is known explicitly: 

. ^an-V^v 7 ^ x > 

X* 0*0 = { A , M pr x n (4-12) 

-4tan" 1 e M v ~ x x <0 
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which consists of two parts of a Sine-Gordon soliton. We remark that this soliton is 
permitted only because the discontinuity allows the field to vanish at infinity. Otherwise 
there could be an infinity of other classical solutions. 

Inserting this into eq. (|4.7|), it gives the estimate for the Wilson loop: 



W, 



c 



N 
1 

N 



32tt 2 



d 6 x 



N 



2N 



( X ' a - vT - M (N - l)(cos( X a/2) - 1) 



-a a A 



(4.13) 



with string tensions a a of 



a. 



g 2 M N-1 j2 
32Ntt 2 2 V N 



8tt N ViV 



(4.14) 



This shows that confinement of quarks exists in this theory for generic values of the 
coupling constants and for finite N, while to extract the behavior with large N it is 
necessary to perform the average for M and then for the Wilson loop j^zr J2 a e~~ aaA . 



The N dependence of M, that is of £, is not known explicitly. As in [14]], £ is the one 
loop partition function in a single monopole background, before the integration of the 
zero mode translation coordinate: 



AT9/2 



rn 



3/2 

w 1 



q 



A(X/g 2 



_ N UlWplc(\/Ng 2 ) 



(4.15) 



The condition of validity of the saddle point approximation, which represents the low 
Debye density, is £ << 1. It can be seen to hold for finite N from the above expression, 
where the exponential vanishes asymptotically. In the limit of iV large one needs a more 
precise estimate on the whole average £. This we will attempt in the next sections. 



4.5 Determinant 

In the previous sections, £(q) plays the part of the statistical-quantum weight of a par- 
ticular background configuration, so that the higher it is, the higher is the importance of 
that particular configuration, although it may have large action. 

Up to now we assumed £ to be some fixed quantity. Now, in order to draw some 
conclusion about the string tension a, we need to find something more precise on it. 

The evaluation of £ is the problem of calculating the functional determinant of the 
fluctuating fields around the one monopole solution. It would be a hopeless problem to 
calculate it exactly in an arbitrary external field, as it is equivalent to the solution of a 
Schroedinger or Dirac equation in an external potential*. 

We will try to extract some information from the high N analysis of the problem. 

*In the expression has calculated numerically with the heat kernel method for SU(2), but the 
estimate of the behavior with N has yet to be done. 
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The idea comes from the fact that after gauge fixing, and better in the unitary gauge, 
the Higgs field has only N components, while its effective action, upon integration of the 
gauge sector, is of order iV 2 . Hence the saddle point should be applicable, and the Higgs 
field is a semiclassical quantity with respect to l/N, which acts like h — > to suppress 
its fluctuations. 

We will treat the fields in one loop approximation around the one monopole configu- 
ration Ap, 4> adding the fluctuating fields a M , tp, so that: A^ = + a M , = + tp. 

In doing so, we are faced with the problem of gauge fixing, because there are zero 
modes of the action. The gauge invariance involves the total field (A^ + a^, <p + tp), and 
one can split the gauge variation between the fluctuating and the background fields in an 
arbitrary manner. 

Among the possible (infinite) choices, one can assign the whole field variation either 
to the fluctuations or to the background. The latter choice is of little or no utility, the 
first, instead, is quite convenient in that it keeps away the gauge invariance problem from 
the background fields. 

So we will keep the background fields in some fixed gauge, and consider the gauge 
group as acting on the sole fluctuations: 



The gauge for the classical fields is left for now unspecified, even the radial gauge 
satisfy automatically the background gauge. In the next section instead we will choose 
for all fields the unitary gauge. 

The partition function in a one-monopole background, £, is a gauge invariant object 
but not gauge independent, at least in principle. Every gauge fixing could provide different 
physical insight, as happens already with spontaneous symmetry breaking. 

According to 't Hooft and Polyakov £ it is better calculated in the so called background, 
"natural", gauge, for two reasons: one is the presence of zero modes, to appear in a 
moment, for the action of fluctuations, and they are best treated in the background 
gauge; the second is that candidates with opposite features, the unitary gauges, usually 
addressed to be non-renormalizable, can be cured only introducing a nonpolynomial term 
in the action. 

This can be seen from the Vandermonde determinant appearing in the measure of tp 
after elimination of the (SU(N)) 'angular' part of tp: Dtp — > Dtp ]J X A 2 (tp(x)). 

There are two ways to deal with this A 2 ; it can be reabsorbed in the measure by a 
(nonlinear but nonsingular) change of tpi variables, but it yields a nonlinear model on a 
singular curved target space of the kind (for SU(2) for example): 



Alternatively, thinking to a lattice, one can exponentiate the determinant and convert 
it in the divergent logarithmic potential: 5(0) I l°g \4>i ~ (fij\d 3 x, but the continuum 



u gauge 

5gauge(a fl ,tp) = S(A„,<t>) (<fy, V) = (^ tt > «]) ■ 



(4.16) 



L = tp-I(d,tp) 2 



(4.17) 
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limit seems problematic due to the ultralocal nature of A 2 . This procedure can be justified 
with care starting from the kind of gauges ['tH,V] |42[] or with other kind of regular 
gauge fixing. 

I can note however, that this divergent extra potential, being exponentiated without 
the help of ghost fields, carries an h factor, and thus already is a quantum correction to 
the action. 

As the theory is renormalizable for every gauge R^, in the effective action there 
should be, then, a divergent term from the massive gauge fields which cancels it. In effect 
it is what we find after the one loop analysis. 

Hence for the lare N saddle point analysis on the effective potential, there should be 
no serious problem from this term. 

We will adopt in the following paragraph the background gauge and in the next the 
unitary one. 



4.5.1 Background Gauge 

Here we take, as gauge fixing: 

F K = - wg[$, tp] = B(x) (4.18) 

In the limit k — > oo the second term reproduces the unitary gauge. This is a variant of 
the pure background gauge (that has B — 0, k — 1), already considered by || or of the 
"natural gauge" that appears in [14]]. It has a long history in the literature, due to its 
double features of renormalizability together with massive ghost fields. 

The usual gaussian averaging of the gauge fixing 5(F k — B), to get a quadratic term 
in the action, is in some contrast with the high N analysis because removes the gauge 
fixing and leaves N 2 degrees of freedom. On the other hand fixing strictly the gauge in 
the case k ^ 1 presents some subtleties. 

We will denote together the fluctuating fields a^(x) and <p(x) with $ = and the 

full quadratic action will be, in matrix form, S qua d r = t $A4.<£>. Scalar product $ • $ will 
imply integration on space-time. 

About the action we just need to say, for now''', that it is annihilated by the following 
zero modes: 



Explicitly it is: 



Squadr — NtX 



([Dp, a,] - [D v , a M ]) 2 + lG av [a a ,a v ] + {g[a u ^ + ip} + if]) 2 + ^V'^ip 



Together with the proper mass term for gauge fields, ^Ng 2 [(f>, a u } 2 , there is also the bilinear mixing 

N([a u ,<p}[D^J} + [avfafivM)- 
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translation = \W = 9 — 

The complete field <3> can thus be decomposed in zero modes plus nonzero ones, $„, 
eigenf unctions of the action: 

$ = R i $® +$ (a) + Z n $ n 

The translational modes (|4.20|) are written in a gauge so that they satisfy the gauge- 
fixing ( |4.18| ) with k = 1, £> = 0. This is useful because they are orthogonal to the gauge 
modes. They are normalized to 




|$«| 2 =Ni = Ntr J + (A0) 2 - (4.21) 

All these zero modes are treated with the standard Faddeev- Popov method to extract 
the integration on collective coordinates Ri, a(x): 

d 3 RVa 5(F k - B) JJ S ($ • . 

W 

Insertion of this unity in the functional integral replaces the flat fluctuations with the 

1 /2 

right variables Ri, a(x) and gives the two determinants. The first Det gives Ili M > the 
second is the Faddeev-Popov determinant 



Det, 



d 
dR~ 



(0 



Det 



5a(x 



Det 



5a(x 



Det[M 



FP\ 



Det 



Kg 



(4.22) 



After elimination of the first delta, the remaining fluctuations represent the functional 
integral restricted the non translational modes: 

Z = e- 5 ^ j d 3 R]jAfl /2 J P<M(F fe ($) - B)Det[M FP }e-f t * M * (4.23) 



The Faddeev-Popov determinant can be evaluated, in one loop approximation, in the sole 
classical background fields. 

Then, as seen from ( |4.21| ), each factor in Ili A/"/^ 2 is the action in a space-time direction 
without the potential. Because we consider spherical monopoles only, all A/i are equal 
and 



,r N 
Mi = ytr 



1 



-G 2 + (D<P) 2 



(4.24) 



They coincide with the action in the BPS limit, hence, after the discussion of section 
PH in the large N limit we also take M = N^-^. 
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One can perform a functional integration over B (with e 2 J- 82 ) to remove the 5(F k 



B) 



Z = J d^e-^A/^Det - ng 2 [^ [0, •]]] / P$e~ S^iM+^M^ 



d 3 i?e- 5 ^ 3 / 2 Det iD^t)^ - m? 2 [0, [0, •]] 



Det 



-1/2 



M + ^M gf 



(4.25) 



Here with *$.Mg/<& = F%(<&) we denote the standard gauge fixing term arising from F K , 
while with D<3> = VaVip, as with the determinant Det, we integrate on translation-fixed 
fluctuations. 

This is analogous to the formula of Polyakov |TJ|], but does not show the explicit 
dependence on k; one has to go through the loop expansion in gauge bosons and ghosts 
in the external field, and proceed to resum all contributions of generic order in g if one 
wants to control the limit k — > 00. Actually in this limit the ghost masses become large so 
that they decouple, but at the same time the coupling with gauge and Higgs also diverges, 
so that the resumming is nontrivial. We will see the first order in an other way in the 
next section. 

Now we proceed in a different direction, we eliminate the 5(F K — B) by direct integra- 
tion of the gauge modes. 

This point of view has the advantage of showing the physical degrees of freedom $ n , 
together with the explicit dependence on k that we want to compare with the unitary 
gauge. 

Briefly, instead of Det 



Det 5/ [M] Det 2 



%M FP 



M. + 2-Mgf we get the decoupled product of two determinants 

plus a measure jacobian dependent on k. The second determi- 
nant cancels formally the FP determinant, and A/2 disappears in the normalization, k of 
course remains in the measure, but the limit is nonsingular 



Let's expand as before the fields as $ = + $o( a ) + We can eliminate the 

zero modes, taking into account the jacobian from $ to (a, Ri). 
We get (still ignoring contributions at more than one loop): 



Z = J d 3 Re~^'l[Afl /2 ^D^t 



DuDp-g* 



]]} jv(^ n )e-^S^ M ^ (4.26) 



The dependence on k and B seems disappeared, but the eigenfunctions $ n have to satisfy 
the gauge fixing so that they are sensitive to k and B. 

Remembering that $ n was satisfying the natural gauge ( Q4.18 ) B = 1, k = 0), to 
pass to a choice with B 7^ 0, k 7^ 1 we have to perform a gauge transformation of the $ n , 
so that they satisfy the new gauge fixing: 



$ n -> $ n + $ (a n ). 



(4.27) 
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This gauge transformation is: 

a n (B,K-l) = \D^D^-g% [0, •] 



(ig(K-l)[^,<p n ]+B) 



(4.28) 



Notice that a n (B, k — 1) — > oo when k — > oo. 

Because we are just mixing with components along the gauge zero modes, the quadratic 
form in the action is unaffected. 

The new basis $ n is still orthogonal to the translational modes. For this reason we 
can retain the same measure of translations A/" 3//2 . 

However we have now eigenf unctions which are non normalized and not even mutually 
orthogonal, and we have to re-normalize the measure. 

m^ip = v[^ n ] =nd£ n -> p[a($„+$oK))] =n d ^- J 



J = \fdet mn [l mn + 0o(ttm)0o(«n)] 



(4.29) 



s _ J(*$+* ( Q ))A1(*+*o(a)) = 9 ~M nn £ 



The new jacobian J carries the information about the gauge dependence on k of the 
functional integral. J, although a formal expression, is explicitly function of the gauge 
parameter. The partition function is: 



f = / d 3 Re- s ^ 2 JDet [b.D, - g% [0, •]]] J JJ d^e-"^ (4.30) 

^ V n 

For example in the case k = 1, 5 ^ we have a n (B,K — 1) = a(B) and J = 
det mn [l mn + \(f>o(a(B))\] It does not depend explicitly on eigenfunctions, but when 
K^lwe have: 



J = Jdet [lmn + $0 («m) " $0 («n)] (4.31) 

a n (S, k - 1) = [£)„£)„ - 3 2 [0, [0, •]]] 1 (ig(K - 1) [0, y? n ] + B) . (4.32) 
We can rewrite J explicitly as a function of k as: 



J 



\ 



det 



lmn + (« - IfD^ 



(m) 



A,A,-^, [0,.]] 



and after rescaling (/c — l) 2 , we can write: 



J = \ /det 



(n) 



A 



det 



lmn+ (k _ 1 1)2 f-D^aj, ^^D,-^^,^,.]]^ ^ 



(n) 



(4.33) 



(4.34) 
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The conclusion that we would like to draw from the present calculation is that in the 
limit k — > oo the singular behavior of the FP determinant in ( 4.25 ) has been canceled by 
the gauge bosons. 

We pass directly to the unitary gauge, then, and we extract some information about 
the large N limit. 

4.5.2 Unitary gauge 

In this section we fix the gauge by requiring the Higgs field to be diagonal. 

Choosing the unitary gauge eliminates all the local gauge invariance apart from the 
little group and the aforementioned discrete Weyl group. 

In the case that we consider, we assume the background Higgs field to have all different 
eigenvalues and the fluctuation to be small with respect to it, so this is precisely the case. 
We will verify a posteriori whether this vacuum configuration is preferred. 

The residual f/ Ar_1 (l) abelian gauge invariance has to be cured in a second moment 
by means of a further gauge fixing. 

In this unitary gauge the Goldstone bosons are explicitly "eaten up" with the standard 
mechanism by the relative gauge fields which acquire one polarization more together with 
the mass. 

At the same time two things happen: first the gauge fixing requires, through proper 
handling of the integration measure, the introduction of a Faddeev-Popov determinant; 
second the massive gauge fields have a Proca propagator, which carries bad behavior at 
large momentum. 

For this last peculiarity the unitary gauge is usually addressed as non renormalizable, 
because the gauge fields produce, even at one loop, a new set of counterterms not present 
in the original lagrangian. 

There is a lot of literature on the self-canceling of some of these non-renormalizable 



divergences, starting from fr4"?|| , the remaining divergences are found (see and therein) 
to vanish on the equations of motions, so that on-shell amplitudes do not suffer of this 
problem. 

We will see that the Faddeev-Popov determinant participates exactly to render the 
theory manifestly renormalizable. 

This can already be inferred from the £ — > oo limit of the i?g gauge (for the charged 



gauge fields): for any value of £, the theory is renormalizable |27|, [43], [44], [45|, [461, so that the 
only counterterms needed are of the same form of the lagrangian. The limit £ — ► oo is well 
defined for the massive gauge propagator, so that assuming some suitable regularization 
(Dimensional regularization is not suitable for this scope because the divergences we have 
are of the kind 5(0) which vanish identically: / d d p = 0) the cubic divergences cancel 
order by order between ghosts and massive gauge fields (or, alternatively, between the 
Faddeev-Popov determinant and the gauge fields). While the charged ghosts acquire an 
infinite mass, thus decoupling, also their coupling to the Higgs becomes large, leaving a 
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correct counterterm to the gauge divergences. It follows that in the careful limit we do 
not expect any new effective interactions. 

A cancelation of this kind has been proven to happen in an abelian gauge model by 
Appelquist and Quinn [fSfl . 

Explicitly we impose the unitary gauge by the constraint^: 

F{<p) = tp char a ed = (4.35) 

It has no derivatives and acts on a field which transforms locally under the gauge group, 
hence it requires the pointwise Faddeev- Popov jacobian: 

5 



-F{S a <p{x)) 



= [<p,-}S(x-y) = [<f ) + V ,-]5(x-y) (4.36) 

F=0 



8a(y) 

that gives the following functional determinant: 

Bet[M FP } = n n a +<Pi- h - = n a2 (^(^)) 

x i<j X 

where A(0) = Oi<j(0i — 4>j)- Here a regularization has to be implicit to make sense of 
the infinite product. 

However Det[.M fp] can be exponentiated without the help of ghost fields thanks to 
the relation DetA = e Trlogj4 to yield an effective potential for the Higgs field: 

Then, because the action is multiplied by a —1/K factor, we need to multiply by a —h 
factor, so that it ends up describing a one loop correction to the bare action (in the form 
of a repulsion of the <j) eigenvalues, as in matrix models). 

In three dimensions this correction has a cubic divergent coupling constant and a non 
polynomial structure. 

Together with this one loop correction, we have then to consider the other contribu- 
tions from the propagating fluctuations, namely the one loop diagrams of gauge fields 
and diagonal Higgs field <p, to calculate, at one loop: 



■fit is the same thing to consider the £ — > oo limit of the background gauges of the last section, 
namely [(f) , 0] = or [0, (p] = 0. 

§Now we would like, later, to exploit the fact that in the unitary gauge the Higgs field has only N 
components, while its effective action is still of order N 2 . In the large N limit this implies that the 
fluctuations tp are suppressed, and <f> is in all respect a classical field. 

In particular it will be possible to apply the saddle point method on its effective action: 

e -m = J Pa ^ e -i[s c! +s, uadr -M(o)/d^io g A 2 W (x))] _ (4 4Q) 

Of course this calculation is still not simple, because it depends on the classical field A. 
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To this aim, we recall the action of fluctuations 



'quadr 



iVtr 



([Ai, a„] - [D„, a^}) 2 + -G^[a^ a u ] + (g[a^ + <p] + y]) 2 + 

+V^"(0Vl (4.41) 



understanding that cp is diagonal. 

All the effect that we want to discuss arises from the massive gauge fields circulating 
in a loop, so we calculate the divergent part of it. 

Because the gauge propagator is constant at large momentum, the loop contributes 
with an arbitrary number of insertions of interaction terms. 

The cubic interactions in the gauge or Higgs fields do not enter at one loop, and 
simultaneous interactions of two a M with a Higgs field plus a background gauge give 
perturbative corrections to what we need, so we leave them apart. 

From the above action ( |4.41| ) we take the relevant quadratic interaction terms of the 
fluctuating ajf with the external fields: 



(4.42) 



where the Higgs field = 4>+<p has been decomposed in a different way: the asymptotic 
constant field 0°° which regulates the gauge bosons mass, plus the remaining classical 
nonuniform background and fluctuating fields which have to be treated as a total external 
field <f r : = 0°° + <f'. 

The loop consists thus of the same charged ajf field run- 
ning along the loop with propagator 



9nu 



p 2 — mfj 



+ 



mfj p 2 



m. 



(4.43) 



it 



g 1 0°° — 0°° | ) and an arbitrary number of inser- 

2/ 




(with mo- 
tions of Vi 

We insert this v%j at zero external momentum, because we are dealing with the diver- 
gent part. We have thus, for n insertions, 



, tJ - <rW - 20f + ^ - - 40- 



Kj 



d 3 ptr 



p2 _ m 2 



+ 



p X p 



mfj p 2 



mf 



*(0)V 



U \'"u. 



(4.44) 



There is a combinatorial factor which comes from the (n — 1)! ways to insert the 
interactions compared with the n! ways to attach the resulting counterterm. Hence the 
factor is l/n. 

Summing up all these divergent contributions we reconstruct the logarithmic potential: 



M(0) En i M^r+^x^)) 



ns(o) (log |0, - 0,i - log |0°° - if i) 



(4.45) 
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It clearly cancels only part of the Faddeev-Popov Vandermonde determinant above in 
( |4.39| ), and leaves the second term function of the Higgs vacuum only: 



-M(0) J logA 2 (0°°) (4.46) 

Now, on one hand it is just a constant which is the same for all the topological sectors 
and can be absorbed in the normalization for what regards £ and all fluctuations, on the 
other it can be thought as a potential correcting the vacuum constant value 0°°. 

This result needs some discussion. A similar mechanism is well known from the study 
of the perturbative corrections to the Higgs potential |42 |: the quantum corrections keep 



a nonzero Higgs v.e.v. also in the limit of no bare breaking [i — > 0. 

From this point of view it provides a repulsive potential for the eigenvalues of the 0°° 
v.e.v., justifying the assumption of (pi ^ 4>j for the vacuum value. 

It may appear strange that we have obtained an effective potential with a divergent 
5(0) ~ A 3 constant, because it seems to be stronger than any of the renormalized other 
terms in the action, in the continuum limit. 

Nonetheless the fact that it depends only on 0°° does not allows us to treat it like 
the other terms in the effective action. The scalar potential V(<j)) is of a very different 
nature because it depends on the fluctuations also. It can lead the field to attain its 
minimum, which has to be the vacuum for the system, in order to be stable against local 
perturbations. 

The effective potential (|4.46|) is expected to play a role let's say, just at the moment 



of symmetry breaking, when the sources should decide which direction in the Cartan 
space to choose. It is at this stage that the Vandermonde potential is important and the 
quantum theory requires one unique direction with nondegenerate Higgs field. 

Eventually one must leave to external sources just the discrete choice among the 
N(N — 1) possible vacuums related by Weyl. 

All this analysis is independent of the large N limit, but for our dilute gas of mono- 
poles, we must would like to know 0°°. In the next paragraph we will find it according 
to the above discussion, for large N. 

4.6 The Higgs vacuum 

The vacuum field plays an important role in the dilute gas picture, because it decides 
if the monopoles are relevant to confinement. 

According to the discussion about the unitary gauge in last section, the vector of 
eigenvalues (f)°° G R^ -1 is defined by the minimum of 

logIC "CI ( 4 - 47 ) 

i<j 
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with the constraint 



oo\2 



(4.48) 



and we recall that p? is of order N, so that the components of 4>°° are of order 1. 

The solution for finite N, although existing, is not illuminating. We instead turn to 
the large N limit and introduce the (non standard) density of eigenvalues: 



N- 



di 



(4.49) 



It is of order 1, as a consequence of last constraint. 

We want to solve for it to obtain the distribution of eigenvalues of the vacuum Higgs 
field. The equation for p(x) is, from (|4.47|), (|4.48|): 



N-P [ a -^-dy = 2Xx (4.50) 
J -a x — y 

where the right hand side is the Lagrange multiplier for the constraint ( |4.48| ) which is 
equivalent to / x 2 p(x)dx = p. 

In the large N limit the Lagrange multiplier is negligible (we do not expect an at- 
traction of eigenvalues coming from this constraint, because p 2 ~ N) and we have to 
solve: 

P f -E^-dy = (4.51) 
J -a x — y 



The by now standard method [5(J is to introduce the resolvent of p as 

F(x) = pj a ^Z~ d y (4-52) 

which has the following properties: It is analytic out of the cut [—a, a] on the real axis; it 
goes to zero at infinity as A; It is real on the real axis [—a, a] excluded; near the cut it 
has zero real part and a discontinuity in the imaginary part given by the unknown np(x). 
The unique function with these requirements is 

F(x) = C J_ (4.53) 
\Jx l — a 1 

from which we finally read (and normalize) the distribution p(x): 

p(x) = \ 2 . (4.54) 

7r vcr — x z 

The result is thus an inverted semicircle law. 
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Its domain is denned by the constraint (|4.48|) (we have introduced the fixed scale 
~p 2 = p 2 /N): 



N x p{x)dx = Nfl 



which gives: 



r.2 



9 

2 ^ ~2 

a = —p 



(4.55) 



(4.56) 



p(x) of ( |4.54| ) represents the following Higgs configuration in the large N limit: 



' 2 -/2cos(vr^). 



(4.57) 



Of course there are N\ equivalent configurations related by Weyl permutations. We show 
the ordered ones in the picture. 



-a 




Eigenvalue distribution p((j)) 



Ordered eigenvalues of (]) 



Let's make a few comments on this result: 

• First, in the large N limit we find that the Higgs eigenvalues remain of order one, 
while in the usual picture there is no indication and one could also assign all the vacuum 
expectation value p to a single eigenvalue. 

• Second, the Weyl degeneracy is completely broken by (f)°°, and no symmetry remains 
apart from the abelian U(1) N ~ 1 . 

• Third, and more important, our vacuum value QL57|) shows a peculiar characteristic, 
namely that, near the edges of the distribution, differences of eigenvalues are of order -j^ . 
As a consequence the masses of gauge bosons in the SU (2) sectors near the two ends are 
vanishing as l/iV 2 .In the middle one has the normally expected 1/N masses. 
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This facts means that there will be an infinite number of very light monopoles, with 
masses vanishing as 1/N. 

The implications on the semiclassical picture of monopole gas are interesting, and we 
will deal with them in the next section. 

On the other hand, independently of the monopole gas, the presence of an infinity 
of massless modes can lead to a new phase of the model, of course present only in the 
N = oo sector. 

Let us mention also a nice correspondence with the solvable cases ofA/" = 2— ►A/"=l 
supersymmetric gauge theories in four dimensions. 

In the SU(N) generalization of the (A/" = 1)SYM±, as analyzed by Douglas and 
Shenker [p2"] , the </>°° field is represented by the points in the moduli space of the J\f = 2 
theory which become the vacua of the theory broken to M = 1. 

The exact solution shows that those points are exactly of the form of our . In 
particular the Weyl degeneracy is completely broken, and moreover the different charged 
sectors will have masses of gauge bosons and of monopoles given again by cos(7r-^) — 
cos(7r-^), for example ~ 1/N 2 near the edges of the distribution. 

We remark that in the supersymmetric theory the ground state appears in a full 
nonperturbative and geometric way, while in our case it must be stressed that we are not 
dealing exactly with quantum corrections, but with a global factor in the Higgs measure, 
expected to be important just during the formation of the system. 

A main feature, the hierarchy of masses, is the same, and this is at least curious. 



4.7 The hierarchy of masses and monopole degeneracies 

Here we come to the observation that in the large N limit there is a hierarchy of masses 
which arises from the model, once given the distribution of Higgs eigenvalues. As A" gets 
large, many of the objects become very light, potentially bringing in a different physics, 
like happens in other cases in this limit. 

We have the gauge bosons and monopole pseudo-masses 

10 • q 



%(q) = #10°° ■ q| S cl (q) = 47rATi 
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which are function of the SU(2) root q which specifies the color sector. Explicitly as a 
function of the Higgs eigenvalues 

. q | = 3-L|#o (4.58) 

when q is the charge in the sector (ij). 

The possible values are N(N — 1) and we see that according to the distribution of 
eigenvalues in the Higgs field, they range in the interval (fi/N 2 , 1). 
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Moreover at the ends of the Higgs domain, ±/x, where the eigenvalues concentrate, we 
have masses of order 1/N 2 , while from the standard differences we get ~ N 2 masses of 
order 1/N and more. 

We can look at the distribution of differences: 

<r(d) = J J dxdyp(x)p{y)5{\x - y\ - d) 

r-d/2 i 
= 4 / dt . (4.59) 
J -« ^{a?-t 2 ){a 2 - (t + d) 2 ) 

which is logarithmic singular for d — > showing the phenomenon: 

aid) ~ — - log a/d. (4.60) 

because masses are proportional to d this is also the distribution of masses in the model 
near the edges. 

The monopole pseudo-mass S c i has a N factor which compensates this vanishing in 
the middle part of the distribution, but still many monopoles tend to have zero action in 
the limit. 

4.8 Dilute gas in the large N limit 

We will use these informations in the dilute gas ensemble of section |4.2| . 

We have to perform the average of £ (q) in the subspace of the simple roots which are 
not orthogonal to a given weight m a . These are of the form: q = m a ± mp {(3 ^ a). 

Then 

t(m*) = -j^-— £ £(m a - mfl) (4.61) 
For large N the average £ has the form: 



9/2 

h 

9 

which leads to the following average 



9/2 

£ ~ N 9 / 2 / dxp(x)^Ae- N ^ (4.62) 

J Q i 



n"*a [ a dx J 1 i(--f) 9/ \-m^ (463) 

J -a \/a 2 — x 2 g 6 

where the large N limit has to be performed. 

Here the first observation is that the exponential in the large N limit becomes a delta 
function (times a 1/N factor). 

But then the (x — b) 9 ^ 2 prefactor drives the integral to zero. In fact the integral is 
expanded and factors of 1/N appear leaving derivatives of the delta function. 
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However we know that also the factor A depends on N, being currently an open 
problem. 

In case A will turn out not to grow faster than iV 2 , the string tension will be suppressed 
and vanish in the large N limit. This could be explained noticing that also in the large N 
limit, as in the BPS case, all the Higgs fields become massless. This important peculiarity 
distinguishes between the large N YMH and pure KM theory at A > 0, whereas for SU (2) 
the Higgs field just provides an explicit scale for the monopole stability. A complete 
discussion will be available by the time the nontrivial quantity A in the one monopole 
determinant will be studied. 



A A remark on the Gribov problem 

Here we want to make some remarks on the Gribov ambiguity in the present case where 
the unitary gauge is chosen. 

It turns out that no Gribov problem is present, and the unitary gauge is thus just a 
feasible gauge with some unusual poperties. Moreover the absence of unphysical fields, 
and the fact that the Higgs field is already diagonal, makes it a good tool to investigate 
the quantum theories. 

The main point observation in the Gribov problem is that the Faddeev-Popov deter- 
minant carries the information on how good is our choice of gauge. It represents the local 
jacobian for the change of variables from the connection space to the quotient by the 
gauge group action where we have our physical theory defined. 

At points in functional space where the FP determinant vanishes, it means that we 
are not taking a complete quotient, in fiber bundle language the section we have chosen 
becomes tangent to the fiber along some gauge direction. 

These points usually mark what is called the Gribov horizon, the name because it is 
the boundary of the maximal region where a section can be continued. 

Here, to be concrete, we have in three dimensions a gauge connection and a matter 
field so the space of fields T is that of couples (A^, <fi), functions from R 3 to the gauge 
algebra. As we have seen T is disconnected in components according to the total magnetic 
charge, as described in section |2|. 

In every component acts the gauge group Q = {g(x) : R 3 — > SU(N)}, and no bound- 
ary conditions have to be imposed at infinity because by homotopic arguments gauge 
transformations connected to the identity do not change component. 

This is best and more appropriately seen in the regular gauges, while in the unitary 
gauge we know that A^ becomes singular along some Dirac string (still with S c i < oo). 
Nevertheless the unitary gauge has nice properties, namely the FP determinant and the 
gauge fixing only depend on the Higgs field, so that one can draw some conclusion. 



The unitary gauge introduced in section [4.5.2|, has a Faddeev-Popov determinant 
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which turns out to be 



X ij 

Now we recall, as remarked at the end of in section |2|, that in nontrivial sectors of the 
Higgs winding at infinity, has necessarily some coinciding eigenvalues at some point. 

This is proven in any regular gauge but is valid also in other gauges because it's a 
gauge invariant statement. 

In the semiclassical picture of monopole gas, at each monopole location two (or more) 
eigenvalues coincide. 

Hence we find that in the unitary gauge the FP determinant above seems to vanish 
identically for any nontrivial configuration. More suggestively one can think that the 
Gribov horizon is made of monopole configurations. The same happens at each field 
configuration throughout the whole nontrivial sectors. 

This would mean that the unitary gauge is an ill defined section, in that it does not fix 
the gauge. Explicitly it would leaves intact some gauge subgroup at the points in space 
where we have a monopole. 

However there are two points which solve this seemingly bad problem. 

• The determinant in the form above is ill defined. It requires us to live in a distribution 
space, whereas we usually consider smooth functions. After this remark it appears evident 
that no gauge invariance remains unfixed, because the smoothness constrains the gauge 
variation at the "origin" to follow that in the neighbor. So there is no such thing as the 
gauge variation at a single point, even if the theory has local invariance. 

• In subsection [4.5.2| we proved that the Faddeev- Popov determinant and its vanishing 



is canceled by the gauge loops, so that if the jacobian vanishes it is just because the change 
of variables is singular, and in fact at the same time the integral takes care of this and 
diverges by the same amount so to correct the measure. 

This is different from the standard Gribov phenomenon where the vanishing of the 
FP determinant is a unavoidable problem in the context of a renormalizable theory. 
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